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> . Abstract 

X: 

■ We calculate the \jvr? c corrections in the inclusive semileptonic widths of 

D mesons. We show that these are due to the novel penguin type operators 
that appear at this level in the transition operator. Taking into account the 
nonperturbative corrections leads to the predicted value of the semileptonic 
width significantly lower than the experimental value. The worsen the 

situation or at the very least, within uncertainty, give small contribution. We 
indicate possible ways out. It seems most probable that violations of duality 
are noticeable in the energy range characteristic to the inclusive decays in the 
charm family. Theoretically these deviations are related to divergence of the 
high-order terms in the power expansion in the inverse heavy quark mass. 



1 Introduction 



Recently a QCD based approach to calculation of total inclusive decay rates of 
heavy (B and D) mesons was developed |j] - |J. The approach is based on the 
systematic expansion in the inverse heavy quark mass within the operator product 
expansion (OPE) 0. In this note we will discuss application of this formalism to 
the calculation of the total semileptonic inclusive width of D° mesons. Unlike the 
previous model calculations the OPE-based method gives us full control over all 
relevant parameters in theoretical expressions. 

The leading perturbative (0(a s )) || 3 and nonperturbative {0{l/m 2 c )) || cor- 
rections have been found previously. It turns out that under a reasonable choice 
of the c-quark mass the predicted value of T S \(D) is significantly lower than the 
corresponding experimental number |37J]. We calculate the next-to-leading nonper- 
turbative correction of order l/m? c and show that it only worsens the situation, or at 
the very least, gives small contribution within uncertainty limits. We then indicate 
possible ways out. 

At the level of C?(l/mj?) terms there arise penguin diagrams generating new, 
four-fermion operators of dimension 6. The penguin graphs were introduced 20 
years ago [|ll|] in the strange particle decays where it was crucial that they produce 
right-handed quarks. In the D meson semileptonic decays the origin of penguins 
is quite different - they appear at the level of the transition operator and give rise 
to contribution of the annihilation type. Usually it is believed that the latter is 
suppressed by chirality arguments. The suppression is lifted, however, due to the 
fact that penguins produce the right-handed quarks, much in the same way as in 



Ref. ITT 



Since the \jw? c terms do not eliminate the discrepancy between the theoretical 
prediction for T S \(D) and experiment a natural question immediately comes to one's 
mind: what went wrong? In estimating the Z)-meson matrix elements of the four- 
fermion operators we use factorization. In the limit of large number of colors, N c — > 
oo, this approximation becomes exact. One may suspect, however, that at iV c = 3 
deviations from factorization are substantial. Can these deviations be a solution of 
the problem? 

Although logically this possibility is not ruled out a priori it is hard to believe 
that this is the case. Indeed, if the problem is to be solved in this way not only the 
matrix element of dimension-6 operators must be enhanced by a factor of ~ 3, its 
sign has to be reversed as compared to what one obtains within factorization. 

The second logical possibility - an enhanced contribution coming from dimension- 
7 operators - also seems very unlikely. 

Thus, we are inclined to conclude that the failure of the standard rriQ 1 expansion 
in the case of r s i(-D) is due to the fact that the charmed quark mass is too light 
for duality to set in. This assertion will be explained in more detail in Sect. 4. 
Here we only note that duality is one of the crucial elements of the calculation of 
the inclusive widths within the heavy quark expansion. Theoretically the onset of 
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duality is related to the behavior of high-order terms in the I/tuq expansion, the 
divergence of the 1/tuq series. The fact that the OPE-based power expansions are 
actually asymptotic is well established |T2[. Very little is known, however, about 
specific details of the divergence. 

An indirect although a very strong argument that the charmed quark is only at 
the border, or even below the boundary, of the duality domain comes from consid- 
eration of the lifetime hierarchy in the charmed family (for a recent discussion see 
[p"3|). Although 0{m~ 2 ) and 0(m~ 3 ) terms qualitatively reproduce the observed 
pattern some of the predicted lifetime ratios (which span an order of magnitude!) 
are off by a factor ~ 2. The predicted (D(m~ 2 , m~ 3 ) deviations from the asymptotic 
limit are typically smaller than what is observed experimentally. Needless to say 
that asymptotically all lifetimes are equal. 

As was mentioned above, the issue of the inclusive semileptonic D decays was 
addressed in the recent literature more than once. The approach to the problem 
accepted e.g. in Ref. is inverted. It is assumed that the theoretical prediction 



for T S \(D) truncated at the leading order of perturbation theory and at the leading 
order of the l/m c expansion (i.e. keeping only 1/m 2 ) is accurate enough to use it 
to fit the values of the quark mass and other theoretical parameters from r s i(Z^) cxp . 
The value of the charmed quark mass emerging in this way is unrealistic. At the 
same time the average value of the heavy quark kinetic energy /i 2 remains essentially 
undetermined. We, instead, use the best available scientific estimates of m c and fi 2 . 
We will see that the results are in direct disagreement with the experimental data. 

Organization of the paper is as follows. In section 2, we discuss the current 
situation. The naive parton result is augmented by its perturbative - to a s - and 
non-perturbative - to 1/m 2 - corrections. Section 3 shows the situation with \jm\ 
non- perturbative corrections. We will see that they don't improve the match with 
experiment. Finally, we discuss ways out of the dilemma. 



2 The starting point 



The theoretical expression for the inclusive semileptonic width of the D meson (c — > 
slv transition), including the leading perturbative and non-perturbative corrections, 
has the form 



192^' cs| 



1 + A {l) a s - ^ - 

2m 2 2m 2 



(1) 



Here m c is the charmed quark mass, and we neglected the strange quark mass; V cs is 
the corresponding CKM matrix element. The coefficient of the 0(a s ) term has 
been known for many years, see Ref. [14]] whose authors merely adapted the QED 
radiative correction to fi —>■ evv (the original QED calculations are published in Ref. 



15|). The explicit expression for A^> depends on what definition of the quark mass 

euu implies the use of the 



m c is accepted. The straightforward borrowing from /i 
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so called pole mass. Although this parameter is not well-defined in full QCD (see 



Ref. |TJ| and the discussion below) it is admissible for a limited technical purpose 



of presenting the 0(a s ) correction. Then 

A W = - — (<,1*-^). (2) 

The leading non-perturbative correction in the l/m c expansion for T(D — > luX s ) 
was first calculated in Ref. M . The term of the first order in 1 / m c is absent , the so 
called CGG/BUV theorem ||, |J. At the level 0{m~ 2 ) the correction is determined 
by two parameters, 

(m^cn d) g 

2M D 4 V D DJ ' 

2 _ {D\^Df\D) ^ 
*** ~ 2M D ' [ } 

Before proceeding to numerical estimates it is worth discussing the parameters 
in Eq. (|I]) in more detail. Numerically the most important parameter is the quark 
mass since it enters in the fifth power. As was shown in Refs. |§ - || it is the current 
quark masses, not the constituent one and not the hadron mass, that appear in the 
systematic 1/rriQ expansions. If we limit ourselves to the 0(a s ) expression for T and 
do not ask any questions about higher-order terms we are free, of course, to express 
the result in terms of the pole mass or in terms of the running mass normalized at 
any point we like - this will merely redefine the coefficient in front of a s in a certain 
way. We are aimed, however, at better accuracy; in particular, we want to include 
the power non-perturbative terms in analysis. 

The use of the pole mass is generically inconsistent; only the running mass can 
appear in any OPE-based expression, see [pT6| . The question is what normalization 



point is relevant. The full expression for the decay rate, including all terms in the 
a s expansion, is certainly independent of the choice of the normalization point /i, an 
auxiliary parameter in the operator product expansion. For a truncated series - and 
we are forced, of course, to truncate the series at the level of the leading, or, at best, 
the next-to-leading term - the choice of [i is important since under a "natural" choice 
the coefficients in the neglected part of the series are small while under "unnatural" 
choices they can be abnormally large. In Ref. |T(| it was explicitly shown that the 
natural choice for m c is /i = Const. m c (see also [[17], [18]]). The leading operator in 
the expansion in the problem at hand is c(i p) 5 c. By adopting the normalization 
point fi = Const. m c we avoid any large corrections. Non-perturbative effects enter 
through the matrix element of this operator; they are also represented by matrix 
elements of other (subleading) operators, for instance, c{i p) 3 iaGc. 

Using the equations of motion one reduces the leading operator to m^cc, where 
both m c and cc are taken at fi = m c . We then evolve cc down to a low normalization 
point, fi m c ; the net effect of this evolution is reflected in a factor of the type 
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c(/i, m c ) = 1 + ai(m c / {j,)a s (m c ) + a 2 (m c //i)a^(m c ) + .... which is, anyway, included 
in the perturbative calculation, see Eq. ([[]) []• Once the operator cc is evolved down 
to a low normalization point we useP] the relation || 

cc = cj c + - — -ciaGc ^cvr 2 c + Oillra 1 ) + total derivatives. (4) 

4m^ 2m~ 

The numerical value of the (one loop) pole mass of the charmed quark was 

determined long ago from the charmonium sum rules ||19|| , rnP° le « 1.35 GeV (see 



also 01). A recent advent of HQET |^T] allows one to conduct a consistency check 



of this estimate. Indeed, let us observe that 

m b -m c = M B -M D + /4(^- - -^) + 0{l/m 3 c , . (5) 

where 

- M B + 3 Mb* 
iW B 

and the same for D. Next, for the pole mass of the b quark a very precise evaluation 

m b = 4.83 ± 0.03 GeV (6) 



is obtained in the recent analysis of the QCD sum rules for the T system To 
be on a safe side we multiplied the original error bars by a factor of 4. It is worth 
noting that it is very difficult - practically impossible - to go outside the indicated 
limits. The central value of m b above implies m c ^1.33 GeV (provided we accept 
the estimate of Ref. |23| for /x^, see below). The most generous error bars in m b 



and fi* (see below) Plus or are translated in ±70 MeV uncertainty in m c . It seems 
perfectly safe to say that m p ° le lies between 1.25 and 1.40 GeV - one can not imagine 
that the one-loop pole c quark mass is less than 1.25 or larger than 1.40. The HQET 
result m p ° le ~ 1.33 GeV matches very well the QCD sum rule number quoted above. 

principle, this factor must contain also terms of order (/j,/m c ) n due to the exclusion of the 
domain below fi from the perturbative calculation. It is important that the power n starts from 
n = 2, and this term conspires with and 

2 The 0(m~ 2 ) terms in Eq.(4) were derived in Ref. 4. The fact that the operators of dimension 
6 are absent in this expansion can be easily established by using the equation of motion, 

1-7Q 1 . 
— ^ — c = ^ — 

2 2m c 

Then 

_l-7o -1 - 7o l-7o 1 

c — - — c = c — — c = - — -cirire 

2 2 2 Am 2 c 

implies that 

1 i 
5(1 - 7o)c = -^c(tt 2 + -oG)c. 

Moreover, using the equations of motion, again we see that cttqC actually reduces to an operator 
of dimension 7 and we arrive at Eq. (4). 
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As was mentioned, in the purely theoretical aspect, it is more consistent to use the 
running mass in the expression for the total semileptonic decay rate. One may choose 
the so called Euclidean mass ^tj) or the running MS mass evaluated at p 2 = 
—m 2 . Both are close to each other numerically, and are smaller numerically than 
m poie s j nce they are deprived of a part of the (perturbative) gluon cloud compared 

to mr fe , 



m P° le = m ™ d 



2 In 2a, 
1 + - + 

7T 



and m p ° le = mf s (m c ) 



4a s 
1 + ^ + -. 
3n 



(7) 



For instance, a fit in the charmonium sum rules yields Jl9| ml pa 1.25 GeV. Since 



our task is limited - in the numerical evaluation we will not go beyond the first 
order in a s - essentially it does not matter which expression for the width is used: 
the one written in terms of the pole mass or in terms of any other mass from Eq. 

(0)- 

Let us first ignore the correction terms in Eq. ([!]) altogether. Then the naive 
parton-model expression 

OV s) lg27r 3l cs\ 

with m c = 1.4 GeV yields T(D — > luX s ) = 1.1 x 10~ 13 GeV, to be compared with 

T(D -> lvX s ) cxp = 1.06 x 1(T 13 GeV. (8) 

The value of m c = 1.4 GeV is at the upper boundary of what we believe is allowed 
for the pole charmed quark mass. We will consistently push the estimates of T(D — ► 
lvX s ) to the high side and, hence, use this value for orientation. 

If, instead, m c = 1.25 GeV is substituted then the parton-model formula gives 
0.6 of T(D -> luX s ) exp . 

The key point, emphasized in the introduction, is that all known corrections, 
perturbative and non-perturbative are negative. Consider first the 0(a s ) correction 
in Eq. (|]). To evaluate it numerically one needs to know the normalization point 
of a s . Usually it is taken to be \i ~ m c . Again this is the question of higher-order 
terms. Recently it was argued within the Brodsky-Lepage-Mackenzie hypothesis 
25] that the 0(a 2 ) terms are negative and large, so that actually the normalization 
point of a s constitutes a relatively small fraction of m c . Since this effect works in the 
direction of reducing T(D — > ^X,) theor and we agreed to push the estimate to the 
high side, we will ignore the 0(a 2 ) terms remembering, however, that our estimate 
will then lie higher than the actual theoretical prediction. 

For consistency we use the value Aqcd ~ 150 MeV relevant to the one-loop ap- 
proximation. Then a s (m c ) pa 0.30, which is consistent with the recent very precise 
evaluation ||22|| , as well as with the previous analyses of data on deep inelastic scatter- 
ing p4"| . (The corresponding two-loop value Aq CD ~ 250 MeV). Notice that bigger 
values of Aqcd sometimes advocated in the literature only worsen the disagreement 
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between theory and experiment. With this value of a s (m c ) we find that the 0(a s ) 
correction in Eq. ([[]) is equal to —0.24, i.e. it further reduces T(D — > /^X s ) thoor by 
a quarter! 

Let us discuss now the non-perturbative 0{m~ 2 ) terms. The value of \j? g is 
known phenomeno logically, see Eq. (^), 

/!% » 0.41 GeV 2 , (9) 



(small effects due to the anomalous dimension of the chromomagnetic operator |26 
are neglected). As for fi 2 , at present this parameter is not measured, although it is 
measurable in principle [27|. At least two independent lower bounds are established 
P^| , |2"T| , P5| which turn out to be close numerically, 

l4 > 0.4 GeV 2 . 



Moreover, the value of /i^ was evaluated in the QCD sum rules p3| , with the result 

l4 = 0.6 ±0.1 GeV 2 . (10) 
Now we are finally able to estimate the 0(m~ 2 ) terms in T(D — > lvX s ) t Y iem . Again, 



c 

,2 . 



trying to increase T theor we use the lower value of /z 2 ; we then conclude that the 
chromomagnetic and kinetic term contribute to the brackets in Eq. ([!]) —0.31 and 
—0.13, respectively. So collecting everything, we have (to 0(a s ) and 0{m~ 2 )) 

T{D -> luX s ) them = r [l - 0.24 - 0.31 - 0.13]. 

We have less than half of the experimental width! Can the 0(m~ 3 ) terms fix the 
situation? As we will see in the next section the answer to this question is negative. 



3 The 1 /m c corrections to the semileptonic width 

Since the 1 / m 2 terms do not solve the problem of the total semileptonic width, it is 
natural to consider the corrections due to l/m? c terms. They can be calculated in 
the standard way within the heavy quark expansion |T[ - . Below the basic points 
of derivation are sketched. 

We start from the weak Lagrangian describing the semileptonic decays 

£(Ai) = ^V sc O. (11) 

Here 

= {sT lx c){vT»e) , (12) 

and = 7^(1 + 75). Eqs. (0), (|12l ) present the lagrangian relevant to the c — > sev 
transition. 
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Next, we construct the transition operator T(c — > X — > c), 

T = ij d 4 xT{£(x)£(0)} = ]T CiOi (13) 

i 

describing the diagonal amplitude with the heavy quark c in the initial and final 
states (with identical momenta). The transition operator T is built by means of 
OPE as an expansion in local operators Oi. The lowest-dimension operator in 
T(c — > X — > c) is cc, and the complete perturbative prediction corresponds to 
the perturbative calculation of the coefficient of this operator. In calculating the 
coefficient of cc we treat the light quarks in X as hard and neglect the soft modes. 
Say, we ignore the fact that in a part of the phase space the s quark line is soft 
and can not be treated perturbatively. Likewise, we ignore interaction with the 
soft gluons. The presence of the soft quark-gluon "medium" is reflected in higher- 
dimensional operators. 

Once the expansion (|13|) is built we average T over the D meson to obtain the 
lifetime, 

1 

T = Im < D T D > . 

M D 1 1 

At this stage the non-perturbative large-distance dynamics enter through the matrix 
elements of the operators of dimension 5 and higher. There are no operators of 
dimension 4 0j. The operators of dimension 5 has been already discussed. The only 
new operators relevant at the level of dimension 6 are the four-fermion operators of 
the type 

Oq = cTqqTc 

where q generically denotes the light quark field and T stands for a combination of 
7 and t a matrices. What particular combination is relevant will be seen from what 
follows. There are two distinct sources of \jvp? c corrections in the total semileptonic 
width: operators of dimension six arising in the expansion for T, and l/m c correc- 
tions in the D-meson matrix elements of the operator caGc and cc. Let us discuss 
these terms starting from dimension 6 operators in the expansion for T. We shall 
see that since T is a Lorentz scalar quantity, only Lorentz scalar operators can con- 
tribute. Thus, the only operator showing up at this level will be of the four-fermion 
type, (9 6 . 



3.1 The four-fermion operators at order a s and in LLA 

A four-fermion operator appears in T in the zeroeth order in a s from the diagram of 



Fig. 1. The corresponding result can be read off from Eq. (17c) of Ref. | 3~0fl , where 



one must put C + = C_ = 1 and eliminate the color factor of 3 from the numerator, 
Im f <°> = f^\V cs \ 2 {(c.T^Ck - (2/3)c a ^ 5 c k ){s k T tM s i )} . (14) 

07T 
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The expression in the braces includes the left-handed s quarks only. If we use the 
standard factorization procedure (i.e. saturation by the vacuum state) for estimat- 
ing the matrix element of this operator over the D-meson state we get zero for two 
reasons. First, in the factorization approximation corresponds to the annihi- 
lation contribution, cs —>■ lu, which shows up only in D s . Second, even for D s the 
chiral structure is "wrong", and after factorization the matrix element (D\T^°'\D) 
vanishes. 



Therefore, the four-fermion operator appearing in Eq. (|14]) by itself is not inter- 
esting. Let us recall, however, that the normalization point of this operator in Eq. 
(|14|) is = m c , and before estimating its matrix elements we must evolve it down to 
a low normalization point. (Of course, it is desirable to go to \i of order of the typ- 
ical off-shellness of quarks inside mesons. Clearly we can not do this since then the 
perturbative calculation of the coefficient functions becomes meaningless. We will 
make a compromise and evolve down to fi ~ 0.5 GeV assuming that, on one hand, 
the coefficient functions are still calculable and, on the other hand, the factorization 
procedure can be used for obtaining the Z)-meson matrix element.) It is straightfor- 
ward to take this evolution into account in the leading logarithmic approximation 
(LLA); as a matter of fact, we just parallel the standard penguin analysis (see Fig. 
2a). What is crucial is that this evolution brings in new four-fermion operators, of a 
different flavor and chiral structure, whose matrix elements over D° and D + do not 
vanish within factorization. Calculating the diagram of Fig. 2a with the logarithmic 
accuracy we get [] 

im f w = -^|K S | 2 {-^1^ (~ayv + Ic-rVc) , (is) 
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where 

f M = 7 M (! - 7s) 

and t a are generators of 577(3), 

(Here A a are the Gell-Mann matrices.) Notice that the light- quark current, 97^°? 
is actually ^D U G^ U ; it includes all light-quark flavors and both, left-handed and 
right-handed fields. This is a typical feature of the penguin contribution [11], leading 
to a non- vanishing contribution of Eq. ( |T5"|) to the D-meson matrix elements within 
factorization. As a matter of fact, we can omit the left-handed part of the light- 
quark current, since, as was explained above, after factorization the term with the 
left-handed part of the current will vanish, i.e., 

E W^Q = E \ {^^q + qf^q} ^ W ■ (16) 

q q Z A q 



3 This result can also be extracted from Eq. (20) of Ref. BC 
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3.2 Full O(aS) calculation 



Sj 

Unfortunately, lnm 2 / /i 2 is not a very large numerical parameter and, hence, neglect- 
ing non- logarithmic terms may seem unjustified. Therefore, instead of summing up 
all logarithms in LLA (which can be readily done, though), it seems reasonable 
to limit oneself to the 0(a s ) calculation including both the logarithmic and non- 
logarithmic terms, the more so that we need the result only for the purpose of 
orientation. We want to convince ourselves that the contribution of dimension-6 
operators to T(D — > X s lv) is negative. 

It is most convenient to carry out the full 0(a s ) calculation using the background 
field technique. There are two versions of this technique - the first one was exploited 



in the context of the inclusive semileptonic decays e.g. in Ref. [31 1, the second 



version, based on the Fock-Schwinger gauge, is reviewed in Ref. [jJ2]. Both versions 



can be applied in the case at hand. Here, the latter one is more suitable. 

The dimension 6 operators come from two sources. To see this, we write explicitly 
the expression for the transition operator. 



f = J d"xe-^Q{x)T^S q {x^)T u L, v Q^) 

Where L^ v is the lepton loop, and S q stands for the light quark Green function. 

As mentioned, in the expression we have the propagator for the s quark. When 
expanded in the Fock-Schwinger gauge this propagator is 

-im 2 K\{m\J — x 2 ) m 2 x 

47T 2 \/-X 2 4:7l 2 X 2 



+ 2 -g^- 2 {2K {m^{-x^m))D a G apl(3 - (D a G a ^ + 

mKi(m\/(— x 2 )) 
x^x a D 7 G a ^ - 3ix^x a D^G a/3lpl5 ) V = + .... (17) 

where m is the s quark mass (needed for infrared regularization). We will see that the 
first of the 1/mj? corrections will come from inserting the second (free) term of Eq. 
(|17|) into T. The second correction comes from terms of order DG in the propagator. 
(Note that the term proportional to G and DG in the s-quark propagator yields 
zero due to their Lorentz structure.) 

The third source for corrections, the matrix elements of dimension 3 and 5, will 
be discussed later. 

We start with the expression for T with the free term of the s-quark inserted. 
This term can be read off from the diagram of Fig. 3. Here we use the expression 
with s-quark mass put to zero from the very beginning since the result is non- 
logarithmic (infrared- stable). Doing the arithmetic, it is easy to show that, in the 
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Fock-Schwinger gauge 



. _ Gl\V cs \ 2 Q{^PQ{0) 
m 384tt 3 

where 

= iD^ - gAp, 

and, 

■A-n — 2 . Q\ X f^P^ g . ^ x p x a(D a Gp fJj ) + .... 

Here we are interested in \jm\ corrections and thus look only at the term pro- 
portional to DG in . Our strategy is to pull A^ to the left since A(0) = 0. When 
doing this, we create commutators which are easily evaluated using the explicit ex- 
pression for Ap. Performing this procedure, we obtain the following contribution to 
the transition operator: 

T|/ree = ~i&F I Ks | 2 ^mfc pt* Oifr pi* Q (18) 

where the subscript 'free' means the free term of the s-quark propagator. Note, here 
we have used the equation of motion 

D a G a atx = -gq^t a q. (19) 
After factorization (see Section 3.4) this yields 

AT = 8a s nf 2 D M D 
T 9ml 

where fo is the axial constant of the D meson. This key constant, fo, is not 
measured accurately enough so far, although some experimental results do exist. It 
seems reasonable to rely on theoretical calculations which were done both on the 
lattice and in QCD sum rules (see Refs. [32] and [33], respectively), thus we choose 
fo — 170 MeV so as to push our estimate for T s i(D) to the high side. We also take 
a s = 0.31, and m c = 1.4 GeV . Plugging in these numericals the above expression 
gives, 

Ar 

— = 0.016. 

Next, we consider the diagram of Fig. 2b. Its contribution to the transition 
operator is calculated in the Appendix. Fig. 2b singles out the DG terms in the 
background field expansion of the quark Green's function, see Eq. (30). The infrared 
cut-off in the logarithm is achieved by ascribing a mass of \x to the s quark line. The 
lepton part of the diagram, which is trivial, must also be inserted, of course. After 
the Fourier transformation we get 
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|dg = iG 2 F \V c 



2 IT/ |2 

cs 



f Q s 2/i m l , 2 \] 

(72^ m «( 1 V + 3 ) | >< 

(2cr a tv + cf, J «v)X 9 - 1 - a r< ( . (20) 

Note that the coefficient in front of the logarithm matches the one in Eq. (15), 
as it should, which was obtained through the logarithmic mixing. Again, after 
factorization, (See section 3.4), we get 

Ar 167ra S/1 m 2 2 w2 „, 

which is, numerically, 

Ar 

= -0.08 

r 

where we have ascribed a value of 0.5 GeV to /i, and used the same values for the 
other parameters as above. 



3.3 0(-M terms from matrix elements of dimension 3 and 5 
operators 



As was mentioned above, the matrix elements of operators of dimension 3 and 5 also 
give rise to the l/m z c corrections in the total semileptonic width. Let us first consider 
Eq. (4). The matrix element of the first term is exactly unity. The matrix elements 
of the second and third terms do contain l/m c corrections, which we discuss here. 

The point is that Eq. (3) expressing < D\c^aGc\D > in terms of the D* D mass 
splitting is valid only to the leading (zero) order in l/m c . Let us observe that the 
spin splitting yielding M|,* — M|, is determined by the following terms in the heavy 
quark Hamiltonian: 

AH = ^—BB + -^BE x 7T (21) 

where B and E are the chromoelectric and chromomagnetic fields, respectively, 
B = gB a t a and E = gE a t a . To the leading order, 

A D = -{Ml. - M 2 D ) = -<ffB>= 0.405 GeV 2 

where < • • • > by definition means {2Mo)~ l < D\c ■ ■ ■ c\D > . At the level of 
l/m c the second term in the heavy quark hamiltonian becomes important in Ad, as 
well as the second order iteration in (2m c ) _1 < oB >. Assuming that both effects 
are of the same order of magnitude, we can roughly estimate the matrix element 
(2m^ 1 ) < BE x 7T > as the difference between A D and A B , 

| (2m c )~ 1 < BE x 7T >|< A D - A B w 0.04 GeV 2 (22) 
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Next, observe that 

1 1 1 

-ca^G^c = -caBc caE x ttc ctDiEAc (23) 

2 m c 2m c 

The last term in the above equation reduces to the four fermion operator which we 
can take into account explicitly. The second term will be estimated as an uncertainty 
in the expression relating < |crG > to A D , 

H% =< % -oG >= A D ± 2(A B - Ad) - (2m c )- 1 4vra s < c llx t a cq ltl t a q > . 

Using factorization for the Oq term above, and the same values for the parameters 
as above, we get +0.01 for the contribution of Oq so that 

H% = A D ± 2(A B - A D ) = 0.42 ± 0.08 GeV 2 

As for /i 2 , we will assume that the error bars in Eq. (10) give the estimate of the 
l/m c part in /i 2 . Moreover, it was shown previously that the sign of the l/m c 
correction in /i 2 is negative (see Eqs. (41) and (42) in Ref. [26]). 

3.4 Factorization and estimate of the matrix element 

We must estimate the matrix elements of the transition operator over the D-meson 
state. As was mentioned above, to this end we use the factorization procedure. We 
realize, of course, that it is not exact - deviations from factorization definitely exist 
- still it seems safe to say that it gives a reasonable estimate of the four-fermion 
operators, especially as far as the signs are concerned. 

The relevant operators are first rearranged, through the Fierz transformations, 
into the form cY^qqY^c, with appropriate coefficients. Then the matrix elements of 
the latter operators are found by saturating with the vacuum intermediate state, 

< D\cT a qqT a c\D >= f D M 2 D , (24) 

< D\cT a t a qqT a t a c\D >= (25) 

Here we used the definition 

< 0\qT a c\D >= if D p a . 

In Eq. ( p5[ ) we accounted for the fact that the matrix element of the color current 
between the vacuum and D meson is zero. 

The L — R chiral structure in Eq. (20) is crucial. We face here an exact analogy 
with the usual penguins. Indeed, if we use factorization while estimating the relevant 
four-quark matrix elements, we see that the contribution of the left-handed light- 
quark current is zero, so we can use Eq. (16). The corresponding result is that of 
the factorized transition operator of Eq. (20). 
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Let us parenthetically note that the penguins we obtain here have no relation to 
the penguin operators of the type 

(ct a s) L J2(Qft a q f )L + R 
f 

contributing to charm nonleptonic decays. Although our penguins look similar, their 
origin is completely different from the usual ones. 

After summing up all new effects, i.e. the effects coming from both the expansion 
of T and the uncertainty of Gq and O n we get, pushing things to the high side, 

Ar 

— = -0.06 ±0.06 + 0.03 (26) 

here, the first number is due to the four-quark terms in the transition operator, the 
second is due to the uncertainty of Gq, and the third is due to G n . 

According to our estimates, the uncertainty in Gq and G n is enough to possibly 
make the total contribution of 0(l/m? c ) roughly zero. The sign of the Gq correc- 
tions is undetermined, however, and a minus sign gives a result which worsens the 
agreement with experiment. 



4 Ways out 

We saw in the previous section that contrary to all the hopes, the contribution 
to the inclusive semileptonic width, however exotic it is, does not solve the problem 
of the deficit of the semileptonic inclusive width. There are several possibilities 
which might explain why the general heavy quark expansion fails to reproduce the 
experimental width. 

First, the factorization that we used while estimating matrix elements can be 
suspected. However, the corrections to factorization can be estimated using the 
method of ref. and they seem small. 

Second, a possibility exists that operators of dimension 7 are important. In 
principle, this may happen since the expansion parameter is ~ y^u^/m c ~ 0.7 and 
is of order unity. However, since the correction due to the dimension-6 operators is 
roughly only 10% it seems unlikely that the dimension-7 contribution will dominate. 

At the moment it seems most probable to us that the discrepancy demonstrated 
above is explained by the fact that the family of charm lies below the duality domain. 
The violations of the quark-hadron duality can be viewed as a cumulative effect of 
all high-dimension operators, taken together. Let us elucidate this assertion in more 
detail. 

Constructing the transition operator as an expansion in m~ l we rely on the Wil- 
son operator product expansion. OPE is well-formulated in the Euclidean domain 
where all field fluctuations can be classified as short- distance and large distance. 
Even in the Euclidean domain the divergence of the non-perturbative series in 1/ m c 
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in high orders produces exponential terms of the type exp(— rnj.) which are not seen 
to any finite order in the expansion fL2 |. To get a rough idea of these terms one has 



to invoke instantons or similar model considerations. From the QCD sum rules it 
is known |2(J that these terms are essentially unimportant in the Euclidean domain 
till surprisingly low off-shellness. 

Kinematics of the problem at hand is essentially Minkowskean since we have to 
take the imaginary part of the transition operator at the very end. One justifies an 
OPE-based procedure by keeping in mind an analytical continuation. In the problem 
of the semileptonic width this may be a continuation in the momentum of the lepton 
pair || - one considers the transition operator at such momenta that one is actually 
off the cuts corresponding to production of the hadronic states, in the Euclidean 
domain. The prediction on the cuts is made by invoking dispersion relations, in full 
analogy with what is usually done in the problem of the total hadronic cross section 
in the e + e~ annihilation. In general, one can analytically continue in some auxiliary 
momenta which has nothing to do with any of the physical momenta. This becomes 
the only option, say, in the problem of inclusive nonleptonic widths. 

Whatever analytic continuation is done, strictly speaking the prediction for each 
given term in 1/ m c expansion refers to the Euclidean domain and is translated to the 
Minkowski domain only in the sense of averaging which occurs automatically through 
the dispersion relations. If the integrand is smooth, however, we can forget about 
the averaging, because in this case smearing is not needed. This is what happens, 
in particular, with the total hadronic cross section in the e + e~ annihilation at high 
energies - the quark-hadron duality sets in and the OPE-based consideration yields 
the value of the cross section at a given energy, locally (without smearing). At what 
energy release is the integrand smooth and can the terms in the l/m c expansion can 
be predicted locally? The existing theory gives no answer to this question. It may 
well happen that at a given (Minkowskean) energy E the deviations from duality 
fall off only as a power of 1/E. Such a regime takes place, for instance, in a model 
discussed in ref. |j36| . This model is definitely relevant for the large N c limit. It 
seems more likely, however, that in the real QCD, with N c = 3, the violations of 
duality fall off exponentially, exp(— (E/Eq) 1 ), and the rate of this fall off is correlated 
with the divergence of the high-order terms in the power series |T2 . 



Very little is known about this aspect of QCD at present, the issue definitely de- 
serves further study which is clearly beyond the scope of the present paper devoted 
to an applied problem. The genuinely theoretical approach would require deter- 
mining the rate of the divergence of the high-order terms in the power expansion. 
Alas, we can not do that, and the 1/mj? term found gives no hint on this divergence 
whatsoever. 

In the absence of theoretical considerations we are forced to rely on phenomeno- 
logical information. The total hadronic width in r decays is a problem close in 
essence to that considered here. Moreover, the r mass is only slightly higher than 
m c . A detailed QCD-based analysis of r decays has been carried out in ref. pTfl , and 
the discussion of the results of ref. in the context the issue of duality is given 
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is ||36|| . There are good reasons to believe that the deviation from duality in the r 
hadronic width is at the level of 7%. Is it then reasonable to expect that descending 
from m T to m c we get a deviation at the level of factor of 1.5 or 2? 

The fact that the onset of duality is not universal, generally speaking, and de- 
pends on the channel considered, is known for a long time [J38|| . Where specifically 
the difference lies between r — > uX and D — > lvX s (apart from the obvious differ- 
ence, m c /m T ~ 0.7) remains to be found. 

Finally, let us emphasize that all attempts to determine parameters of QCD or 
HQET from the analysis of the heavy quark expansion in the charm family must be 
viewed with extreme caution and are hardly reliable in view of the uncontrollable 
theoretical situation discussed above. 
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5 Appendix 

In this Appendix we discuss the calculation of the transition operator associated 
with Fig. 2b. We use the Fock-Schwinger gauge x^A^ = and the background field 
method (see Ref. |32| for details). We shall need the part of the propagator S(x, 0), 
see Eq.(17), for the quark with mass m in external gluon field that contains odd 
number of 7 matrices and is proportional to DG. Note that in the limit m — > the 
propagator becomes singular, 



AS -> \n(m^(-x 2 )) (27) 
The contribution of diagram of Fig. 2b to the transition operator is equal to 

T Fig . 2b = Im|Vg 2 ^ J e^cr M A^(x,0)r i ,c(0)(tr[ 1 S(0,x)r^(x,0)r^) (28) 

where the trace term in the brackets represents the lepton loop, and AS* stands for 
the DG part of the s-quark Green function. Here S(x, 0) is the free massless quark 
propagator: 

S(x, 0) = x/(2ttV) (29) 

The direct calculation uses the formulae of Ref. [39| for the imaginary part of the 
Fourier transform of the integrals of the type 



K n (m^J(-x 2 ))/(-x 2 ) p d 4 xe ipx . (30) 
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We just write Eq. fl2~T| ) explicitly in x space using the propagators from Eqs. (|T 
and (|29|) and then use the formulae of Ref. |39| to convert the resulting expression 



into the imaginary part of the integral in Eq. (PS[). We then see that the diagram 
of Fig. 2b gives the following contribution to Im T: 

~3^^V ^ + 2pV)c(0) 7 ,(l + 75 )£> a GV(0) (ln(^/A 2 ) - 5/6) (31) 

We now go the reference frame connected with the center of mass of the heavy quark 
where p u = p 13 = m c (l, 0, 0, 0). Note also that for heavy c-quark 

c l0 D a G a0 c = c(0) lp D a G af3 c (32) 

The Eq. Q3"2|) follows from the fact that cy^c ~ 0(l/m c ) We finally obtain the 



following 'piece' of the transition operator 

Im T Fig . 2b = G 2 F \V cs \ 2 {^m 2 c (\n(m 2 c /A 2 ) - 5/6) 

(33) 

2(cTpt a c + c 1 pt a c)q lp t a q 

The logarithmic term comes from the singularity, Eq. (p7|). This is not all however. 
Part of the contribution to the transition operator of Eq. (^) comes from the 
infrared domain, and its contribution has nothing to do with the contribution we 
are interested in. The infrared contribution is given by the contracted loop of Fig. 
4. It is easily calculated, giving contribution 

Im f Fig . 4 = G 2 F \V cs \ 2 {^m 2 c (\^ 2 /A 2 ) - 3/2) 

(34) 

(cTpt a c + c~if3t a c)q^pt a q 

Subtracting Im T Fig 4 from Im T Fig 2 b we get the transition operator Im T DG of Eq. 
(20). 

Figure Captions 

Fig. 1. The four-fermion term in the transition operator as it appears at the 
level oP s . 

Fig. 2. (a) The penguin graph for the four-fermion operator in T. (b) The 
diagram with the DG term in the s quark line. 

Fig. 3. The diagram with the free s quark line determining Tf ree . 

Fig. 4- Subtraction of the infrared part from the coefficient of the operator 
cDGc, see Fig. 2b. 
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